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Abstract. Using the HST archive WFPC2 observations and rotation curves measured by Vogt et al. (1996), we 
constructed self-consistent light and mass distribution models for three disk galaxies at redshifts z = 0.15,0.90 
and 0.99. The models consist of three components: the bulge, the disk and the dark matter. Spatial density 
distribution parameters for the components were calculated. After applying k-corrections, mass-to-light ratios for 
galactic disks within the maximum disk assumption are M/Lb ~ 4.4, 1.2 and 1.2, respectively. Corresponding 
central densities of dark matter halos within a truncated isothermal model are 0.0092, 0.028 and 0.015 in units 
Mq/pc^. The light distribution of the galaxies in outer parts is steeper than a simple exponential disk. 
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1. Introduction 

The study of dark matter halo properties and central den- 
sities allows to constrain possible galaxy formation mod- 
els and large scale stucture formation scenarios (Primack 
2002; Khairul Alam et al. 2002). It is interesting to anal- 



yse the corresponding data at different redshifts. For this 
kind of analysis it is needed to know both the distribution 
of visible and dark matter. Unfortunately, the structure 
and mass distribution of stellar populations, and there- 
fore, the discrimination between visible and dark matter 
in very distant galaxies is still largely unknown. 

General morphological evolution of galaxies at inter- 
mediate and high redshifts has been studied, for example 
by Sm all et al. (|1997| ), Andreon et al. ( |1997| ), Lubi n et a l. 
( |1998D , Fabricant et al. (^000[ ), van den Bergh et al. (|2000D , 
van Dokkum et al. (2001a), La Barbera et al. (2002), Im 
et al. ( |2002 ). Several statistical correlations and color gra- 
dients in clusters at redshifts z = 0.3 — 1 have been stud- 
ied by Stanford et al ( |199^ ), Ziegler et al . ( |1999| ), van 
Dokkum et al. ( |2000| ), SagUa et al. ( |2000| ), Tamura & 



Ohta (2000). More sophisticated correlation parameters 
for elliptical galaxies, the Fundamental Plane parameters 
involving dynamical information, have been determined 
by Kelson et al. (|1997D , van Dok kum e t al. ( [1998| , |2001b| ). 
Kelson et al. (|2000|), Treu et al. (|200lD. 



These studies indicate that most of the stellar content 
of early- type galaxies has formed at redshifts z > 2—3, and 
has passed normal chemical evolution and there is no sig- 
nificant age difference between field and cluster galaxies. 
Within the standard ACDM model, field ellipticals must 
be younger than cluster ellipticals. On the other hand, the 
fraction of early-type galaxies in clusters is smaller and the 
fraction of E-l-A (post-starburst) galaxies is higher at in- 
termediate redshifts, indicating some morphological trans- 
formations from z ~ 1 to the present time. Theoretical 
models indicate that truncation of star formation rate can 
successfully explain this kind of morphological changes 
(van Dokkum & Franx, |200l| ; Bicker et al. |2002| ). 

Less is known about the structure of spiral galaxies at 
intermediate redshifts. For analysing the emission line pro- 
files of disk galaxies, it is needed to distinguish the rota- 
tion and dispersion components. Spectrograph slit width 
is usually only slightly smaller than galactic dimensions 
which makes derived rotation curves uncertain and often 
only maximum values of rotation velocities can be esti- 
mated. 

By now the maximum values of rotation velocities 
have been measured for distant galaxies up to redshift 
1.34 (Colless et al. 
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1994; Forbes et al. 199e 



1997; Simard & Pritchet 199? , Schade et al 
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1999 



& Stanford |2001|). When plotted to the TuUy-Fisher di- 
agram, these measurements allow to compare the mea- 
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sured galaxies with the local ones. Rotation curves of very 
distant galaxies have been measured only by Vogt et al. 
(|l99|, \m% , Ziegler et al. ( pOO^ ) and Steidel et al. ( pOO^ . 
These data allow more detailed modeling of the structure 
of corresponding galaxies, especially if the curve extends 
beyond the solid-body rotation part. 

To construct a self-consistent photometrical and dy- 
namical model of a galaxy, both surface photometry and 
rotation curve data are needed. In the case of galaxies 
at intermediate redshifts the typical scale is 5 — 10 kpc 
per arcsec and in order to determine the parameters of 
galactic components, it is recommendable to use also high- 
resolution HST photometry in addition to other photo- 
metrical data. 

For our modeling, the galaxies had to match the fol- 
lowing criteria. Their velocity profiles had to be undis- 
turbed, symmetric with respect to the galactic centre and 
with sufficient extent. Galaxies seen almost face-on were 
not suitable due to uncertainties of the inclination an- 
gle of the galactic plane. We also could not use galaxies 
being seen edge-on, because for these galaxies the sur- 
face brightness distribution is largely influenced by galac- 
tic absorption. For detailed modeling, we selected three 
galaxies with rotation curves measured by Vogt et al. 
( |1996D : 074-2237, 064-4412, 094-2210. For deriving the sur- 
face brightness distribution of these galaxies, we used im- 
ages from the HST archive. General properties of these 
galaxies are given in Table |^. In the present work we take 
Hq — 65 km/s/Mpc and go = 0.5. 

2. Photometrical structure 

The images of all three galaxies were retrieved from the 
archive of the Hubble Space Telescope (WFPC2, ob- 
servational programs 5090 and 5109, filters F814W and 
F606W). For processing we used pipeline-reduced on- 
the-flight calibrated images. 12 exposures in both colors 
were available for the galaxy 074-2237 and 4 exposures 
for the other two. Details about image processing and 
derivation of surf ace b rihtness profiles are described in 
Tamm & Tenjes (2001). Here we summarize only basic 
steps, bringing out especially these aspects which differ 
from the previous processing. Processing was done with 
IRAF/STSDAS packages. 

In order to elliminate cosmic rays, exposures were com- 
bined using the STSDAS task combine (option crrej). 
Before estimating the background, combined images were 
filtered to smooth down noises (task adaptive). We found 
no evidence of galaxy distortion caused by adaptive fil- 
tering. As the investigated galaxies were rather small in 
CCD frames, it was sufficient to subtract in all cases a 
constant background value from the images. Background 
values were determined by interpolating the values of the 
surrounding regions of the galaxies. This stage differs from 
the one wc have used earlier. New on-thc-fiight calibration 
of the HST archive data has improved the quality of the 
original images, especially in the case of the galaxy 064- 
4412. The background levels were (in WFPC2 intensity 



counts) 11.7 (I)and 19.8 (V) for the galaxy 074-2262; 6.40 
(I) and 7.25 (V) for the galaxy 064-4412; 6.25 (I) and 7.25 
(V) for the galaxy 094-2210. Images of the three galaxies 
after background subtraction are given in Fig. |l| 

The resulting images were deconvolved using the Lucy- 
Richardson algorithm for modeling PSF with the help of 
TinyTim 6.0. V-band images of 064-4412 and 094-2210 
were not deconvolved, because due to a low signal-to- 
noise ratio deconvolution process gave unrealistic results. 
In order to calculate isophotes, task ellipse was used. The 
fitting process was carried out several times for each of 
the galaxies in order to estimate and reduce uncertainties. 
The scatter of the dots in Fig. |^ illustrates uncertainties 
of isophote fitting. 

Before final calibration, we estimated the corrections 
of the luminosity flux for the geometric distortion and the 
charg e trans fer efficiency effect of the WFPC2 (Holtzman 
et al. 1995a). We added 2% to the signal to compensate 
for the geometric distortion of 064-4412 and 094-2210 and 
1.5% in the case of 074-2237. To transfer the luminosity 
counts into standard Johnson magnitudes, the calibration 
described by Holtzman et al. ( |l995b| ) (H95) was used. For 
the filter F814W, we used the fiight photometric system 
Eq. (8) and the coefficients from Table 7 in H95: 

/ = -2.5 log (DN/t) - 0.062(1/ - /) + 0.025(V" - lf + 

20.839 -I- 2.5 logiGR) + 0.1 + 0.05. 

For F606W filter, we used the synthetic calibration Eq. (9) 
and coefficients from Table 10 in H95: 

V ^ -2.5 log {DN/t) + 0.254(1/ - /) + 0.012(1/ - 1)"^ + 

22.093 + 2.5 log{GR) + 0.1 + 0.05. 

Here GR is gain ratio as defined in H95, the term 0.1 
corrects the formula for infinite aperture and 0.05 is the 
correction due to long exposures. V—I was calculated iter- 
atively. Luminosity fiux was corrected for the cosmological 
distortion by a factor (1 -|- z)^. 

Luminosities have been corrected for a bsorption in our 
Galaxy according to Schlegel et al. (1998). 

We estimated the k-correction in two ways. First, the 
observed colors were transformed to such standard col- 
ors for which the k-correction would be minimal. For this 
we used the method described in van Dokkum & Franx 
(|1996D and Kelson et al. (|2000|). At redshift z = 0.15 a 



suitable transformation would be from V to B and for 
redshifts 0.9 and 0.99 from I to B. For all three galaxies 
we used synthetic spectra of a Sbc-Scd galaxy derived by 
Coleman et al. (1980). Zero points relating AB magnitudes 
and Johnson magnitudes were taken from Frei & Gunn 
(1994). Second, the k-corrections were estimated by us- 



ing the chemical evolution models by Schulz et al. (2002). 
Again, at z = 0.15 V color was transformed to B color, 
assuming the rest-frame (B — V) = 0.63, at redshifts 0.90 
and 0.99 I color was transformed to B color, assuming the 
rest-frame (B — I) = 1.76. For the nearest galaxy, we took 
solar metallicity, for the two others Z = 0.008. The results 
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Name" 


RA 


DEC 


Hubble 


2" 


Scale 


Inclin." 


M% 




(2000) 


(2000) 


type 




(kpc/") 


(deg) 




074-2237 


14^17^45 48 


52°27'59.9" 


Sc 


0.154 


2.76 


80 


-20.0 


064-4412 


14'^ 17™ 53.4= 


52°29'41.1" 


Sc 


0.988 


6.71 


68 


-20.6 


094-2210 


14''17™31.3= 


52°26'08.8" 


Sbc 


0.900 


6.65 


60 


-21.5 



^ Galactic names and redshifts are from Vogt et al. (1996), inclinations are from Simard et al. (2002) 
^ Absolute B-magnitudes in galactic rest-frame are from our models (Sect. 4) 




Fig. 1. Images of the galaxies (from left to right) 074-2237, 064-4412, and 094-2210 in F814W color after background 
subtraction. 
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Fig. 2. Surface brightness distribution of the galaxies (from left to right) 074-2237, 064-4412, and 094-2210 in I and 
V colors. Both original and deconvolved profiles are given. 



obtained by two methods nearly coincide in the case of the 
nearest galaxy, for more distant galaxies the agreement is 
within ±0.15™. Averaging the results of the two methods 
yields the following relations: 

B{z = 0.15) = T^ + 0.32, 
B(z = 0.90) =7 + 0.90, 
B(z = 0.99) ==7+1.02. 



Here V , I and {V — I) are observed Johnson magnitudes 
and B is the Johnson magnitude in galactic rest-frame. 
The final galactic surface brightness profiles are presented 
in Fig. ^ by open circles. Error bars characterize uncer- 
tainties of both background determination and isophote 
construction. 
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Fig. 3. Surface brightness distribution of the galaxies 074-2237, 064-4412, and 094-2210 (from left to right) in rest- 
frame B color (open circles). Solid line - surface brightness distribution of the best fit model, dashed lines - surface 
brightness distribution of model components. 



3. Galactic models 

Although several stellar populations can be discriminated 
in spiral galaxies, in the case of modeling very distant 
galaxies it is reasonable to limit the main stellar compo- 
nents to two - the bulge and the disk. Parameters of any 
other stellar component (central nucleus, metal poor halo, 
thick disk) could not be determined at present from avail- 
able photometric observations. To construct a dynamical 
model, a dark matter component - the dark halo must be 
added to visible components. Of course, the photometrical 
and dynamical models must be mutually consistent. 

When handling only surface photometry data, it is con- 
venient to describe the galactic structure with the help of 
some simple parametric formulae for the surface density. 
Most common formulae are the de Vaucouleurs formula, 
the exponential law or the more general Sersic formula. 
Most flexible is the Sersic formula with an additional pa- 
rameter, allowing to con trol t he steepness of the surface 
density decrease (Sersic 1968 ). Density distribution pa- 
rameters must be determined e.g. by the least squares 
method. Using also kinematic data and constructing a dy- 
namical model consistent with the photometry, the same 
density distribution law must be used for rotation curve 
modeling (and for velocity dispersion curve, if possible). 
In the case of the Sersic formula, it can be done for spher- 
ical systems with an integer Sersic index m (Mazure & 
Capelato 2002). For a noninteger index and ellipsoidal sur- 
face density distribution a consistent solution for rotation 
curve calculations is not known. 

In the present paper, the density distribution parame- 
ters are determined by the least squares method and they 
can have any value. Thus, it is not reasonable to use the 
Sersic law for surface density distribution. It is more phys- 
ical to start from the spatial density law. An analytical ex- 
pression for spatial densities, reproducing well the Sersic 
formula, was derived by Trujillo et al. (2002). Although 



for the disk an analyticaly simpler form (|i|), allowing to 
use more easily least-square fitting simultaneously for light 
distribution and rotation curve. 

In our model the visible part of a galaxy is given as a 
superposition of the bulge and the disk. The spatial den- 
sity distribution of every visible component (the bulge and 
the disk) is approximated by an inhomogeneous ellipsoid 
of rotational symmetry with constant axial ratio q and the 
density distribution law (Einasto 1969) 



p(a) =p(0)exp[-(a/(fcao))i/^] 



(1) 



where p{Q) — hM / {'iirqa^) is the central density, M is the 
component mass. Next, a = B? + jcp' ^ where i? and z 
are two cylindrical coordinates, ao is the harmonic mean 
radius 



(2) 



where dM = Airqa^ p{a)da is the mass of the ellipsoidal 
shell of thickness da. Coefficients h and k are normaliz- 
ing parameters, depending on N, which allows to vary 
the density behavior with a. The harmonic mean radius 
characterizes rather well the real extent of a component, 
independently of the parameter N. The definition of nor- 
malizing parameters h and k and their calculation is de- 
scribed in Tenjes et al. ( 1994 ), Appendix B. Eq. fallows 
sufficiently precise numerical integration and has a mini- 
mum number of free parameters. 

The dark matter (DM) distribution is represented by 
a spherical isothermal law 

pia)^l^pmi + &r'-[i + &r') «<«° o) 



a > a 







this formula is useful, we decided to use for the bulge and 



Here a° is the outer cutoff radius of the isothermal sphere. 

The density distributions for the bulge and the disk 
were projected along the line of sight, divided by their 
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mass-to-luminosity ratios / and their sum gives us the 
surface brightness distribution of the model 



L{A) 



2 

1=1 



Pi{a)a da 

hJA W^^W 



(4) 



where A is the major semiaxis of the equidensity ellipse 
of the projected light distribution, Qi are their apparent 
axial ratios = cos^ 7 + 9^ sin^ 7. The angle between the 
plane of a galaxy and the plane of the sky is denoted by 
7. Summation index i designates two visible components, 
the bulge and the disk. 

The masses of the components we determined from the 
rotation law 



Jo 

V\R)^j2vKR), 



Pi{a)o?da 



,2)1/2^ 



(5) 



(6) 



is 



where G is the gravitational constant, e — 
eccentricity, and R is the distance in the equatorial plane 
of the galaxy. Summation is over all three components 
now. 

Model parameters ao, g. Lb, M and N for the bulge 
and the disk were determined by a subsequent least- 
squares approximation process. As the first step, we made 
a crude estimation of the population parameters. The pur- 
pose of this step is only to avoid obviously unphysical pa- 
rameters - relations (Q) and (^) are nonlinear and fitting 
the model to the observations is not a straightforward pro- 
cedure. Next a mathematically correct solution for each 
galaxy was found. Details of the least squares approxima- 
tion and the general mod eling procedure were descr ibed 
by Einasto & Hand (|1989D , Tenjes et al. (|1994 |l99S|) . 

The rotation curves of the galaxies 074-2237, 064-4412, 
094-2210 have a rather small extent (see Figs. We do 
not have any additional information about the dynamics 
of these galaxies at large galactocentric radii. For this rea- 
son, the outer cutoff radius of the DM component remains 
undetermined at present. We fixed aP — 5ao on the basis 
of the structure of nearby galaxies. This value influences 
the rotation curve behavior only in extreme outer regions 
of galaxies and does not influence our results. Then we as- 
sumed by analogy with nearby galaxies that the rotation 
curve of these galaxies remains flat at least up to 30 kpc 
(Sofue & Rubin pOOl| ). 

As we have no information about velocity dispersions 
in the central regions of our galaxies we can calculate only 
the circular velocities (Eqs. If the value of velocity 

dispersions is comparable to the value of rotational ve- 
locity at a particular distance, the latter does not cor- 
respond to the value of circular velocity. The difference 
is known in galactic dynamics as the asymmetric drift. 
Typical emission-line dispersions in galaxies at intermedi- 



model velocities must be higher than the observed rota- 
tional velocities. 



4. Results and discussion 

Calculated from our photometry, the total apparent 
Johnson V luminosity of the galaxy 074-2237 \s V — 
19.28 ± 0.05, which coincides well w ith the Vgoe lumi- 
nosity estimates by Vogt et al. ( 1996 ) 18.70 and Simard 
et al. (2002) 18.78. The mean apparent color index is 
(V — I) = 1.5 ± 0.2. The difference between the Johnson 
V and Veoe is 0.49™, if (V - I) =1.5. 



The calculated total apparent Johnson I luminosity of 
the galaxy 064-4412 is / = 22.43, which also coincides well 
with the /814 luminosity estimates by Vogt et al. (1996) 



22.38 and by Simard et al. ( |2002D 22.33. The mean appar 
ent color index is (V — I) = 2.2. The difference between 
the Johnson I and IgM is 0.12™ if (V - I) =2.3. 

The total apparent Johnson I luminosity of the galaxy 
094-2210 from our photometry is / = 2 1.47. The Ian lu- 
minosity estimate by Vogt et al. (1996) is 21.41 and by 



ate redshifts are 30-100 km/s (Koo et al. 1997; Guzman et 
al. 119971; Im et al. 120011). Thus, within central 0.5" - 0.7" 



Simard et al. (2002) it is 21.21. The mean apparent color 
index is (V - I) = 2.4. 

The final two-component models fit the photometric B- 
profiles with a mean relative error of 0.2-0.5 per cent. For 
the fitting of the rotation curves, we constructed two mod- 
els for each galaxy. One kinematical model corresponds 
to the mathematically best fit at all distances, even in 
the central regions. However, due to the arguments given 
in the last paragraph of the previous section, this model 
is unphysical, because it underestimates the bulge mass. 
Therefore, a second model was constructed with the bulge 
mass fixed at higher values. At present we have no inde- 
pendent basis to estimate the bulge mass. We hope that 
two different models allow the reader to understand how 
the bulge mass variation influences the overall behaviour 
of rotation curves in all three concrete cases. Real bulge 
masses and thus also mass-to-light ratios may be higher 
than the values resulting from Table ||. To derive real bulge 
mass-to-light ratios, additional observations (central stel- 
lar velocity dispersions) are needed. 

The parameters of the final models (the axis ratio q, 
the harmonic mean radius oq, the total mass of the popu- 
lation A/, the structural parameter N , the dimensionless 
normalizing constants h and k, B-luminosities and corre- 
sponding mass-to-light ratios) are given in Table |. Bulge 
and DM masses are given in parentheses to indicate that 
on the basis of present observational data it is not possi- 
ble to determine these values uniquely. The last column in 
Table | gives the masses for low bulge mass models. The 
final models are denoted by thick solid lines in Figs. ||-^ 
Models of individual visible components are represented 
by dashed lines, the model of the dark matter compo- 
nent by a thin solid line. In Table ||, calculated from the 
models, integrated parameters are given: apparent I-band 
magnitude, total mass of the visible matter, total intrinsic 
luminosity and mass-to-light ratio of the visible matter. 
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Table 2. Galactic model parameters. 



Name 


q 


(kpc) 


M 


M/Lb 


N 


h 


k 


Lb 
10'" 




074-2237 




















bulge 


0.7 


1.1 


(0.2) 


2.2 


0.75 


2.459 


0.7822 


0.092 


(0.05) 


disk 


0.1 


7.1 


6.4 


4.4 


0.50 


1.571 


1.128 


1.46 


6.5 


dark matter 




60. 


(170.) 






14.82 


0.1511 




(170.) 


064-4412 




















bulge 


0.7 


1.5 


(0.3) 


0.78 


0.51 


1.592 


1.117 


0.38 


(0.05) 


disk 


0.1 


7.3 


2.7 


1.17 


0.43 


1.397 


1.229 


2.30 


3.0 


dark matter 




40. 


(150.) 






14.82 


0.1511 




(150.) 


094-2210 




















bulge 


0.7 


1.7 


(2.1) 


1.9 


0.67 


2.115 


0.8801 


1.08 


(0.7) 


disk 


0.1 


9.3 


6.4 


1.19 


0.25 


1.080 


1.445 


5.38 


8.0 


dark matter 




60. 


(250.) 






14.82 


0.1511 




(250.) 



" Mass distribution model with low bulge masses (upper panels in Figs 4-6). 
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Fig. 4. Rotation curve of the galaxy 074-2237. Open and 
grey circles - observations by Vogt et al. (1996) at the op- 
posite sides of the galaxy, thick solid line - rotation curve 
from the mass model, dashed lines - rotation velocities due 
to visible components, thin solid line - rotation velocities 
due to the dark matter component. Upper panel (model 
1) corresponds to low bulge mass, lower panel (model 2) 
corresponds to higher bulge mass. 

Table 3. Derived integrated parameters from models. 



Name" 


mi 




Lb 


(m/lu. 




(mag) 


IO'^Mq 


1O'»L0 




074-2237 


17.78 


6.6 


1.55 


4.3 


064-4412 


22.43 


3.0 


2.68 


1.1 


094-2210 


21.47 


8.4 


6.46 


1.3 



200 



> 




r (arcsec) 

Fig. 5. Rotation curve of the galaxy 064-4412. 
Designations are the same as in Fig. |[ 



It is interesting to compare the central surface bright- 
nesses of galactic disks with the results obtained for lo- 
cal galaxies. For the galaxies modeled here, the central 
brightnesses of disks are ^b(O) = 21.4, 21.3, 21.1. The 
mean value for two galaxies at redshifts z = 0.9 — 1 is 
/^s(0) = 21.2. For local galaxies the original value by 
Freeman (1970) was 21.65. However, this value probably 
depends on the galactic morphological type and luminos- 
ity, and thus has considerable scatter (Graham & de Blok 
pool ). For galaxies at redshifts z = 0.5 — 1 Schade et al. 
( 1996| ) derived i^b{0) — 20.0, when transformed to con- 
ventional broadband magnitudes, but also with a rather 



200 



> 



200 - 
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r (kpc) r (kpc) 

10 15 20 5 10 



15 




Fig. 6. Rotation curve of the galaxy 094-2210. 
Designations are the same as in Fig. |[ 



significant scatter. Our values are clearly higher. However, 
it must be taken into account that central brightnesses of 



disks have been determined by Schade et al. (1996) within 
the assumption of exponential disks and r^/"' bulge. When 
fitting the disks of our two higher redshift galaxies with pa- 
rameters = 1, the resulting disk central brightnesses are 
brighter by 0.6 mag, i.e. the mean value is /iB(0) = 20.6. 
Disk central brightness is sensitive also to bulge param- 
eters. Our bulge profile is quite different from the r^/"* 
law. And finally, when taking into account also the un- 
certainties in k-correction, our values for the disk central 
brightnesses are not peculiar. 

Figure]^ shows that the light distribution in outer parts 
of the galaxies studied in the present paper is rather steep 
(the parameter N is small). In all cases, the disk density 
distribution in outer parts is steeper than for a simple ex- 
ponential disk. Central concentration parameters (Trujillo 
et al. 2001) for these galaxies are rather small, about 0.15. 



Steep decrease of the surface brightness in outer parts of 
galaxies may be the result of the outer truncation of the 



disk. According to Pohlcn et al. (200C), the truncation oc- 
curs at distances of about 3 disk scale lengths. In the case 
of galaxies studied here, the surface photometry does not 
extend to so large distances. 

The central density of the dark matter p(0) = 
hM/ {Anag), characterizing the rotation velocity gradient 
at the center caused by the dark halo, is rather well deter- 
mined from observations and nearly independent of the 
choice of the dark halo outer cutoff radius. Calculated 
from our models, the central densities of DM components 
are p{0) = 0.0092, 0.028, and 0.015 in units Mg/pc^. It 



Fig. 7. Rotation curve model of the galaxy 074-2237 with- 
out the dark matter component. Disk masses are 8.2 (disk 
1) and 8.1 (disk 2) in units of 10^°Mq. Designations are 
the same as in Fig. 0. 



is important to stress that these values are calculated for 
maximum disk models and are thus the lower limits. This 
must be taken into account especially in the case of the 
first value (the galaxy 074-2237) where models with lower 
disk mass and higher DM density give equally good fit for 
the rotation curve. The visible and DM masses become 
equal for this galaxy at R — 17 kpc in our maximum disk 
model. A mass distribution model of the ga laxy 074-2237 
was conctructcd also by Fuchs ct al. ( 1996| ) who used for 
mass determination a model of spiral arm formation via 
swing-amplification. They concluded that in order to cre- 
ate a two-armed spiral structure, higher DM densities are 
needed when compared with the maximum disk models. 
According to their model, the central DM density of 074- 
2237 is 0.025. 

It is possible to construct mass distribution models for 
galaxies without assuming the existence of a DM compo- 
nent. Corresponding models are given in Figs 0-^. It is 
seen that the calculated rotation curves start to decrease 
after the last observed point of the rotation curves. Thus, 
within the available at present observations, the models 
without any dark halo fit well with observations. 

As mentioned earlier, due to the limitied extent of ro- 
tational curves it is not possible to determine any param- 
eters of the DM distribution. Figures show that ro- 
tation curves can be well fitted without DM. However, 
we think that there are enough independent arguments in 
support of the existence of DM and thus, realistic models 
must include a DM component. For this reason, on the 
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Fig. 8. Rotation curve model of the galaxy 064-4412 with- 
out the dark matter component. Disk masses are 8.1 (disk 
1) and 7.7 (disk 2) in units of IO^'^Mq. Designations are 
the same as in Fig. 0. 
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Fig. 9. Rotation curve model of the galaxy 094-2210 with- 
out the dark matter component. Disk masses are 12.5 (disk 
1) and 11.4 (disk 2) in units of IO^^Mq. Designations are 
the same as in Fig. H. 



basis of nearby galaxies weumed that rotation curves re- 
main flat up to the distances of about 30 kpc (see Sofue 



& Rubin 2001). This assumption still does not enable to 
determine the total DM mass and its characteristic ra- 
dius. Only the central density of the DM component can 
be found. Therefore, the mass and the radius of the DM 
component given in Table § are simply fixed. It is possible 
to fix, e. g. the DM mass two times smaller, decreasing 
correspondingly the DM radius in a way that the central 
density remain unchanged. The central density of DM can 
be dermined from models within the assumption of the 
flatness of rotation curves. The central density of DM is 
not sensitive to the uncertainties of bulge mass. 

The central density of the dark matter distribution is 
a parameter allowing to constrain possible cosmological 
models and corresponding parameters at early periods of 
structure formation (see eg. Navarro & Steinmetz |2000|, 



Shapiro fc Ilie v |2002| , Khairul Alam et al. |2002| , Zentner & 
Bullock 2002). We think it is reasonable to derive also the 
central DM densities at different redshifts. In the present 
paper we have used only one density distribution model 
for DM (Equation ^). In subsequent papers we intend to 
enlarge the sample of modeled galaxies and study the de- 
pendence of model parameters on different DM distribu- 



tion models (e.g. models by Navarro et al. 1997, Burkert 



1995) 
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